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Calculations of B leptonic decays and B — B mixing using NRQCD heavy and Asqtad light valence quarks on 
the MILC dynamical configurations are described. Smearing has been implemented to substantially reduce the 
statistical errors of the matrix elements needed for the determination of fs- The four-fermion matrix elements 
needed for the determination of /^^ Bb, have been calculated and a preliminary resuh is given. 



1. Introduction 

B — B mixing is a key process in standard 
CKM analysis. Calculations of the B{Bs) lep- 
tonic decay constant fsifB,) and bag parame- 
ter Bb{Bb,) are needed to constrain Vtd{Vts)- 
In particular, we aim to calculate the ratio 
Jb^ \/BbJ JbVBb accurate to the few percent 
level allowing a determination of Vts/Vtd in which 
theoretical errors do not dominate those from ex- 
periment. 

The Asqtad improved staggered formulation 
has allowed precise determinations of a variety 
of quantities T, and recent progress has allowed 
the use of improved staggered valence quarks in 
heavy- light simulations 0- Here we report on 
the progress of B leptonic decay and mixing cal- 
culations using NRQCD heavy quarks and im- 
proved staggered (Asqtad) light quarks. We use 
the MILC Uf =2 + 1 dynamical configurations 

2. B leptonic decays 

In previous work ,4^ the chiral behavior of 
/bsV^^^ib7//b\/^"b was masked by large statis- 
tical errors. Here these errors have been sub- 
stantially reduced with the introduction of smear- 



ing. Two dynamical ensembles were used; one in- 
corporating a light dynamical quark mass rnf = 
fnu,d = fns/4: while the other has to/ = TOs/2 
where to^ is the physical strange quark mass. The 
valence light quark mass niq was varied from 
to nis/S, i.e. some runs were in the partially 
quenched approximation. For b quarks, the stan- 
dard tadpole improved Lattice NRQCD action 
correct through 0(1/ (amt,)^) was used along with 
zeroth and first order in l/arrib currents, a^^ and 
TOf, were fixed by T while rriu^d and nis were fixed 
by TT and K respectively. 

A ground state hydrogenic style wavefuction 
was used to smear the heavy quark at both source 
and sink and the optimal radius was found as that 
which minimized the fit errors while maintain- 
ing a reasonable x^/dof. For each combination of 
source and sink smearing, the functional form of 
the k*^ order in 1/ arrib correlator is given by the 
oscillating function 

G{t)^ J2 C('=^')(-l)^'*e-'"^*. (1) 
j=o 

Bayesian multi-exponential fits were done to the 

(k) 

fc = 0, 1 correlators, and $^ were determined as 
= V2C(fe"). (2) 
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fs correct through 1-loop is given through 



+ (l + as(pi+P2))<U (3) 



where 



(4) 



and the factors po, pi, P2t C have been perturba- 
tively calculated 5 . 

Figure 1 compares values from single 

correlator fits to simultaneous multi-correlation 
matrix (smeared source and sink) and vector 
(smeared source only) fits. It can be seen that 
the errors are significantly reduced in the matrix 
case. 

Figure [21 shows the effects of smearing on 
— ^bJ^b- Note that the most chiral par- 
tially quenched non-smeared result incorporates 
only half the statistics of its smeared counter- 
part, while the statistics match in the other cases. 
Fully unquenched and partially quenched chi- 
ral extrapolations which incorporate staggered 
taste breaking effects 6 are underway to the 
new smeared data, in addition to more fully un- 
quenched simulations at mq/ms = 1/8. 

3. B - B Mixing 

3-point correlators were generated with the cre- 
ation of a 5 at < = —tg, conversion to a i? with 
the use of a four-fermion operator O at t = 
and destruction of the B at t ^ tg. The same 
simulation parameters as in the B leptonic decay 
case were used but only so far with nif — 0.01, 
niq — 0.04 (i.e the Bs). No smearing has yet been 
performed, and so far only zeroth order in l/arrib 
currents have been considered. 

The continuum four fermion operator {Ol)^^^ 
is given in terms of lattice operators at 1-loop by 



where 



= [1 + PLLas\{OL)lat 
+ PLSOls {Os)lat 



(5) 



Ol = [V;q7^(1-75)^,] [^q7m(1-75)V'J 
Os = [^q(1-75)^,] [V5q(1-75)^,] (6) 
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Figure 1. Determinations of a'^^'^^^' using single 
correlator (no smearing), matrix (source and sink 
smearing) and vector (source or sink smearing) 
fits. 



and pLL and pLS have been calculated perturba- 
tively IS). 

The lattice correlators have the functional form 



jtB g-rrijtE 



^ (_l)fc*fle-"''*s. (7) 

A preliminary Bayesian fit has been successful for 
nexp = 4, as can be seen in Figure 13 which com- 
pares the data and fit function of the effective 
amplitude (correlator with ground energy contri- 
butions removed). The oscillations in both ts and 

can clearly be seen to be nicely represented by 
the fit function, and the fit has a good x^/dof. 
Work is underway to obtain fits at other Uexp in 
order to solidify these results. 

{OL,s)iat are determined through 



1 



{O 



A. 



{OL,S) 



3 x-L,s)/ai - — 7^ (8) 

zAlBa-^ k,BB 

where ^bb — C'^°'^'; the 2-point ground state 
leading order in 1/ anib amplitude from the previ- 
ous section. Note that the fit is done directly to 
the 3-point correlators without first taking the ra- 
tio over the 2-point. This means that fits at low t. 
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Figure 2. ^$ vs. ratio of light and strange va- 
lence quark masses. Closed symbols indicate that 
smearing has been used while open symbols indi- 
cate smearing has not been used. Squares repre- 
sent fully unquenched results and triangles repre- 
sent results in the partially quenched approxima- 
tion. 



where the correlator has good statistical noise, are 
possible using multiple exponentials without the 
inefficiency and ambiguity of needing a plateau. 
Bb is defined through 

{0,)^^lflMlBB (9) 

and our results preliminarily give 

Ib, VBsAmb) - 0.f97(f6)(28)GeV. (10) 

The first error arises from the fact that the fits 
are still preliminary and the second error is from 
systematics, the main contributions being from 
the neglection of 0[KQCDl'mb) and 0{a1) con- 
tributions. 

4. Conclusions 

We report here on the status of B leptonic de- 
cays and B — B mixing simulations with NRQCD 
heavy and Asqtad light quarks. The use of smear- 
ing has significantly reduced the statistical errors 



Figure 3. Fit result to four-fermion correlator 
with ground state energy contributions removed. 
Squares show data points while lines show fit 
function. 



for the leptonic decay case. Chiral extrapola- 
tions and more fully unquenched simulations are 
in progress. A successful fit to the BB correlator 
has allowed a preliminary result for fs^ \/ Bb^ ■ 
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